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ABSTRACT 

We discuss the implication of the observation of supernova rehc neutrinos on the study 
of the star formation rate (SFR) in galaxies. The limit recently obtained at Super- 
Kamiokande (SK) is already marginally significant: The SFR we derived ipi^o) < 
0.040MQyr-iMpc~^ (at a 90% CL) is about twice the SFR estimated from radio 
observations, and five times the rate from Ha allowing for uncertainties in the model 
supernova neutrino flux. 
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1 INTRODUCTION 

Neutrinos emitted from stellar core collapse fill the uni- 
verse as a diffuse background radiation. The feasibility 
for the detection of t hese neutrinos has been consid- 
ered by many authors ([ Bisnovatvi- Coga^^^eid0^_^8^ 



Krauss. G lashow fc SchramrnT^S^ 
MalanevI IIQOtI: 




jKa£linghat j_Steigman fc Walker! 



Totani, Sato fc Yoshiif 
Hartmann & Woosley 

booot 

Ando. Sato fc Totani I2002D . While estimates of the 
expected neutrino flux depend much upon authors, the au- 
thors are generally negative regarding the feasibility of their 
detection. The problem is large backgrounds from solar and 
reactor (anti)neutrinos at low energies and atmospheric 
neutrinos at high energies; a possible window in between 
(say, the neutrino energy =15-50 MeV) is masked 
by a large background from electrons produced by decay 
of low energy mu ons that escape d etection in the water 
Cerenkov detector JZhang et alllQQSTl . The decayed electron 
spectrum from muons, however, is precisely known, and the 
Super-Kamiokande (SK) group has demonstrat ed that this 
background can be subtracted llTotsukallioOll) . The limit 
derived on the neutrino event rate in the 18-50 MeV range 
is close to the value indicated by some optimistic estimates 
of the supernova relic neutrino flux, which encourages us to 
scrutinize the problem. 

A major uncertainty in the calculation of the su- 
pernova relic neutrino flux is in the star formation rate 
(SFR) and its evolution towards the past. The work 
with high redshift galaxies over the last flve years, 
however, has provided us with significant insight on 
the evolution of the global SFR. The estima t es include 
the use of UV emissivity jMadau et aJ Il99d iLillv et all 
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. IConnoUv et 


alJ 


1997 


ISteidel et al.lll999llTrever et alJ 


199S 


. ISuUivan et alJ 


200t 


. IWilson et alJ l2002^. Ha fluo- 


rescent emission 


iGaUeeo et al.l Il995l. iTresse et alJ l200l 


Glazebrook et al. 1999. Sullivan ct al. 


2000*), radio emission 



and far-infrared emission ^Flore^taLlTl999^ . Most of the 
estimates of the SFR are convergent to « 0.2 dex among 
different authors, if the same observational techniques are 
used. The large uncertainty, however, resides in which tech- 
niques are to be used; Current estimates show an uncer- 
tainty of a factor of ~ 6 (0.8 dex). In particular, the SFR 
estimated from UV depends largely on exti nction correc- 
tions. Madau et aL. (.199^) took Eb-v = 0.1, ISteidel et all 
(1999) indicated 0.15 and lsullivan et alJ (|2000D derived 0.13. 
This corresponds to an uncertainty of 0.36 dex in the SFR. 
The prime purpose of this paper is to consider whether we 
can obtain any meaningful constraints on the SFR from the 
current observation of the supernova relic neutrinos. 

Another focus of this paper is to derive a lower limit on 
the relic neutrino event rate expected in the SK detector un- 
der reasonable assumptions on the input to the calculation. 
The calculation of the event rate suffers from uncertainties 
besides the SFR. Among the most important uncertainties 
are those in the spectrum of neutrinos. In this paper we 
try to reduce the uncertainty in the neutrino spectrum us- 
ing observation of neutrinos from SN1987A as a constraint. 
We may assume that it is typical of type II supernovae, 
since we expect that the physics of core collapse is similar 
even if the optical appearance may have a wide variety. The 
other uncertainty concerns neutrino oscillation, in which 
and Vfi, partly interchange during propagation through 
vacuum and Earth. Recent neutrino oscillation experiment 
showed unambiguously the presence of neutrino oscillation. 
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and determined the oscillation parameters iFukuda et alJ 
ll998l.lAhmad et alJliooH iFukuda et al.ll200j) . We can now 

calculate accurately the effect of oscillation both in vacuo 
and in Earth. This is no longer a source of uncertainties. 

Throughout this paper we adopt the natural units, c = 
fi = 1, and the Boltzmann constant fcs = 1. 



2 THE LOCAL SUPERNOVA RATE AND THE 
LOCAL STAR FORMATION RATE 

A number of extraealacti c supernova surveys 
Jvan den Berg h fc Mc Clurel Il994l . ICappellaro et al.l Il997l 
iTammann et al.i il994>l yield the local supernova rate in 
units of SNu, i.e., the number of supernovae per 1O^''Lb(0) 
per 100 year for each morphological type of galaxies. 
We translate it into the rate per unit cube of spatial 
volume, by averaging over morphological fractions of 
nearby galaxies E/SO : Sa-Sb : Sbc-Sd = 0.32 : 0.28 : 
l|Fukugita. Hogan fc Feeble jll99?j : hereafter FHP), and by 
multiplying the B band local luminosity densi ty of the 
universe Cb = 2.4 ± 0.4 x 10*/iLq Mpc"^ (iBlanton et aJ 
I2OOII: lYasuda et all 120011^ . where h is the Hubble constant 
in units of 100 km s~^Mpc~^. Counting both type lb and 
Ic in addition to type II as core-collapse supernovae, we 
obtain the supernova rate -Rsn as 1.98 x 10"'', 2.11 x 10"'' 
and 4.44 x 10""' /i^yr"'Mpc"^ for the three surveys. We 
take the geometric mean of the three values and refer to 
the largest and smallest of the three as the allowed range: 
2.65 X 10"* 



10 



similar estimate 4.7 x 10 '^h?yr 



yr 'Mpc This is compared with a 
'Mpc"^ bv iMadau et alJ 



For a given star formation rate tli{t) (in units of M© 
yr"^Mpc"'^) _RsN is calculated as 



RsN = i>{t) 



drruj>{m)/m 
/™" dm(t){m) 



(1) 



where ^(m) is the initial mass function (IMF) of stars, for 
which we take the Salpeter form d>(m) ^ m"''^^ for m > 
IMp, and continue to the IMF of ICould. Bahcall fc Flvnnl 
l|l99ffl for a low mass (see FHP) so that the integral in the 
denominator is extended to zero mass. We set m„ = IOOM0. 
Since the Salpeter IMF has been adopted in virtually all 
literature that derived SFRs, our calculation of SFR can be 
directly compared to those, only with a downward correction 
of a factor of 1.65 that arises from the difference between the 
genuine Salpeter IMF with a lower cut off of O.IMq and our 
prescription. We take the critical mass for type II supernovae 
rric to be between 8 and 10 Mq according to Nomoto ( 1984). 

The local star formation rate derived from Q is 
log ^(fo)(Moyr"'Mpc"^)= -2.09+°-j3, where we take h = 
0.72. This may b e compared with the SFR from an Hq sur- 
vey -2.11±0.04 fTres se et ai]l200l Iclazebrook et allll999l: 
we take their independent data points). 
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Figure 1. Star formation rate (SFR) inferred from Ha(/3) (filled 
squares), UV emissivity (open and filled triangles), far-infrared 
(filled circles) and radio (crosses) observations as a function of 
lookback time. The SFR estimated from the local supernova rate 
is shown by open circle at Gyr. The two dashed lines denote the 
90% SK limits for the case of the minimum event rate (the con- 
servative limit; hatching attached) and for the case of maximum 
event rate within uncertainties of the model supernova neutrino 
flux. The dotted line is fit j^. All data use the modified Salpeter 
IMF described in the text. 



In Figure Q we present estimates for the SFR 
as a function of the lookback time (= to — t with 
to the present time, in units of Gyr). The data are 
taken from Glazebr ook et al. (19 99), Sullivan ct al. (200C|), 
iTresse et a'Lrj2002l) . ISteidel et al.l (1999) and lHaarsma et alJ 
ll2000r ~which cover most of the S FR work to date We include 
in the Figure the UV estimates of lWilson et alJ ^02) (filled 
triangles with thin drawings). These data, however, largely 
disagree with other estimates and the reasons are unknown: 
so we do not refer to them further in this paper. We take the 
cosmology of ilo = 0.3 and A = 0.7 with h = 0.72 to draw 
this figure. The solid points are the SFR from Hq (and H/3) 
with the extinction estimated using the Balmer decrement, 
and the open points refer to the values from UV emissiv- 
ity assuming zero extinction corrections. Radio observations 
are denoted by crosses. The grey points are obtained from 
far-infrared (ISO) observations. The SFR estimate from the 
supernova rate is shown at the zero lookback time. We also 
show the constraint from the supernova relic neutrino obser- 
vations obtained from the present calculation given below. 

The figure shows that the SFR obtained from a single 
indicator exhibits nearly an exponential dependence as a 
function of the lookback time, at least for z < 1, as expected 
in t he closed box model an d also in CDM model calculations 
fe.g. lNaeamine et alj200ll) . For example the SFR for t—to < 
7 Gyr obtained from Ha is fitted well with 



Star formation rates 



log V(t)HQ = (-1.96 ± 0.04) + log h + 0.216/i(to - t). 



(2) 
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The SFR derived from far-infrared observations is consistent 
with this curve. The radio data give a consistent slope, but 
the normahsation is dex times higher. The data from UV 
emissivity also indicate a line parallel to ||5J, but located 
lower by 0.50 dex With th e standard extinction law iSeatonI 
ll97fll.lCardeUi et al.ll989l^ Eb-v = 0.20±0.05 brings the UV 
data consistent with the Ha data. In the following argument 
we assume the exponential law for the SFR 



Table 1. Mean energies (in units of MeV) of supernova neutrinos. 



(3) 



with r = 2.8 Gyr for z < 1. This law may not hold for z > 1, 
but star formation at such high redshift is insignificant in 
a estimate of the supernova relic neutrino flux for Ei, > 
12 MeV, which is of our current concern, and we do not 
need to specify any accurate functional form, as we will see 
below. We take V'(^o) as a parameter. Alternatively, we may 
use Eb-v as a parameter taking logi/'uv(io) ~ —2.45 from 
UV emissivity with zero extinction as a fiducial value. 

Neutrino spectrum from type II supernovae 

We must deal with the uncertainty of the neutrino flux 
emergent from type II supernovae. The dominant neutrino 
emission arises from pair creation in the optically thick ob- 
ject. Hence the luminosity of each species of neutrinos is 
approximately equal, as demonstrated by many neutrino 
transport calculations given in Tabled The total neutrino 
luminosity is close to 3 x 10^^ erg, since the mass of all 
observed neutron stars takes a universal value of IAMq. 
We allot a 20% error to this estimate. The mean energy 
of neutrinos depends on details of calculations, e.g. rang- 
ing from 12 to 20 MeV for Ve- The mean neutrino energies 
satisfy {-Ei^J < {Ei^^) < {E^^^), where v^t includes v^, Vr 
and their antiparticles. Neutrino transport calculations give 
EvfijEi)^ = 1.6; see Table Q We take this ratio to be 4/3, 
but the result of our calculation is not sensitive to this ratio 
once we introduce the observational constraint. 

The model neutrino spectrum from supernovae is some- 
what deviated from the zero-chemical potential Fermi distri- 
bution. It is usually parametrised by introducing an effective 
chemical p otential r\, i.e., / = [exp iE^lT^ -ri) + \\-^ with 
77 = 1 - 3 Jjanka fc Hillebrandt|[l989l) . 

For water Cerenkov neutrino detectors the only reaction 
we must consider is VeP e'^n. The cross section of u^^O 
is > 20 times smaller. There is, however, an important con- 
tribution from i'^ — + Ve due to neutrino oscillation. Recent 
solar neutrino experiments at Sudbury and SK show that 
the mixing between and is nearly maximal. With the 
matter effect the neutrinos emergent from a supernova are 
1^2 — — sin Ovs + cos OVfi and vi — cos + sin Ou^ , which 
are the mass eigenstates. The detected in detectors are 
therefore cos^ 6 times th e Ve flux an d sin^ 6 time s the 
flux where sin^ 261 ~ 0.96 iAhmad et al . 2001, Fuk uda et alJ 
I2OOII) . This in principle increases the neutrino detection rate 
due to higher energies of the flux. 

The matter effect of Earth somewhat modifies the mix- 
ing ratio for the neutrino flux that passes through Earth. 
This effect is calculated assuming Earth as a sphere of a 
constant matter density, pe (Earth) ~ 3.2g/cm^. 

Constraints from SN1987A 

A t the epoch of SN1987A Kamiokande iHirata et all 
119871) and Irvine-Michi gan-Brookhaven Collaboration 
(1MB) teionta et all Il987^ detected neutrinos from core 
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Table 2. 90% confidence allowed range of the neutrino energy 
and temperature from 1MB neutrino events for SN1987A. 



V 


(B^J (MeV) 


T^e (MeV) 





10.5 - 14.9 


3.33-4.72 


1 


11.0 ~ 15.3 


3.31-4.60 


2 


11.4 - 15.8 


3.16-4.39 


3 


12.0 - 16.4 


3.01-4.11 



collapse. The gross characteristics of these neutrino events 
agree with what are expected. Here we use the detection of 
neutrino events at 1MB, which has a larger flducial volume 
and is more sensitive to rare, higher-energy neutrino events, 
to constrain the higher energy spectrum of supernova 
neutrinos. We estimate the event number Mmb at 1MB (5 
kton water) as 



iViMB = 3.3 X 10'^ 



dE^ap{E^)Fi,^{E^)x(E^ 



(4) 



where ap{E^) is the cross section for UeP e^n with the 
neutrino energy E^, x{Ejj) is the trigger efficiency, and Fa^ 
is the flux at the 1MB detector, 



F^^E,) = PlMB(S,.)jri,, + (1 - Pimb(£.))J^^^, 



(5) 



where Pimb{Ev) is the conversion probability for — » t-e 
due to neutrino oscillation including matter effects of Earth 
and J^i is the neutrino flux for species i without oscillation, 
for which we obtain 

E, 



Ti = 3.99 X 10^°cm"^MeV 



r.. \ ( E 
MeVy \^1053erg 



,r/^\ 
VMeVy 

C7(r,)/(r,,£.), 



(6) 



where Exox,v is the total neutrino energy, T^^ is tem- 
perature of neutrino for ith species and C(r\) = 
J dE^E^f{ri,E^)/ J dE^E^f{0,E^). The calculation of 
fi MB is standard and is carried out in a way similar to that 
in lLunardini fc Smirnov| i200ll) . taking neutrino trajectory 
inside Earth for SN1987A. Applying Poisson statistics to 8 
events observed at 1MB, we obtain 90% confidence limits on 
{Ev^) and Tj?^ as shown in Table |5| This constraint removes 
much of the uncertainty of the model neutrino flux: hence 
the result of our calculation in what follows depends orfly 
weakly on the assumptions on parameters we assumed for 
the model neutrino flux emergent from supernovae. 

Supernova relic neutrino flux 

The neutrino flux J^. is calculated as 



f -dt 

dz — {l + z)Ly^ ((1 + z)Ey)R'i^{z) 
dz 



(7) 
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= Etot,,. — 7fiC{ri)f{ri,E^ 



(8) 



where the neutrino luminosity Ci, is given by 
120 

Zf is the formation epoch of galaxies, but the integral is 
dominated by a low redshift region. The supernova rate is 
related with the SFR as, 



7?sn(z) = 1.22 X W'^iP{z)/Mq 



(9) 



for rric = 8Mq. We remark that the flux ||7|l is independent 
of the cosmology, since its dependence in dt/dz is compen- 
sated by the volume factor of tp{z). We use for the ex- 
ponential law for z < 1 and assume a constant for z > 1. 
We take tp{to) as a free parameter. This parameter is also 
translated to Eb-v, taking the SFR from UV emissivity 
with zero extinction as a fiducial. The extinction is written 
A lo gi/i = 2A2Eb-v for z < 1 using th e standard extinction 
law llSeaton|[l979l. ICardeUi et allligsgll . 

The calculation of supernova relic neutrino flux Ja^ (af- 
ter neutrino oscillation effect) is carried out in a similar way, 
with the modification that the Earth effect is integrated over 
all directions. The event rate with a fiducial volume of the 
22.5 kton SK detector is calculated as 



log f (to) 
-3.4 -2.2 -2 -1. 



-1.6 



— , — , — ^ — , — , — , — I — , — , — , — I — , — ^ 
<E-> = 12MeV, <E-.> = 16MeV 




O.S 0,3 
E(B-V) 



RSK 



4.2 X 10" 



VMeV; 



J^.{E^) 



MeV cm^^sec^^ 



,(10) 



where Eu^^.^^^^^^^^ is the minimum (maximum) energy of the 
neutrino detection, which we take i5y„i„ ~ 18 MeV and 
^fmax = 50 MeV. Fig. |5| shows an example result for the 
neutrino event rate for SK as a function of Eb-v for a typ- 
ical parameter set indicated in the figure. It is seen that the 
event rate exceeds the current SK limit, 2.0 yr~^ at 90% 
confidence for the 18 - 50 MeV window llTotsukallioOI') . ff 
Eb-v^OA. From this figure we can read the range of Eb-v 
which is consistent with the SFR from Ha at the zero red- 
shift. The SFRs from Ha and the local supernova rate are 
consistent. 

Fig. |3| shows the relative importance of supernovae at 
different redshifts for neutrino events. The solid histogram 
corresponds to neutrino events with energy between 18 MeV 
and 50 MeV. It shows that half the events arise from low- 
z supernovae {z < 0.25). If we decrease the energy of the 
detector window to 12 < iJ < 18 MeV, non-zero redshift 
supernovae become more important. This means that we 
can learn the evolution of the SFR from gross spectroscopy 
of supernova relic neutrinos. This histogram also shows that 
the contribution from supernovae at 2 > 1 is insignificant 
in so far as our consideration is restricted to iJ > 12 MeV. 
The neutrino spectrum is presented in Figure 2] 

Constraint on the star formation rate 

We derive a constraint on the SFR by requiring that 
the supernova relic neutrino event should not exceed the 
SK limit allowing for uncertainty of the model neutrino flux 
shown in Fig.|5| the result is shown in Figure above. The 
most conservative limit is obtained by taking rUc = WMq 
and minimising the number of events at 1MB for SN1987A. 
At a 90% confidence the limit means 



Figure 2. Neutrino event rate for SK as a function of Eb—v- The 
hatch is the limit from SK. The range between the two dashed 
lines is consistent with SFR from Ha observations, and that be- 
tween the two dotted lines is consistent with the local supernova 
rate. 
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Figure 3. Neutrino events partitioned according to the redshift 
of supernovae. The solid (dotted) histogram corresponds to the 
energy window of 18-50 MeV (12-18MeV). We take the param- 
eter set the same as that in Fig. 1^ 



log ^(fa) < -1.40, or Eb-v<OAS. 



(11) 
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4 CONCLUSIONS 

The limit on the SFR derived from the supernova rehc 
neutrino observation at SK is already marginaUy signifi- 
cant even if we include the uncertainty of the model su- 
pernova neutrino flux. With the current data we can con- 
clude that V(io) < O.O4OM0yr~^Mpc"^, which is 5 times 
the estimate from Ha (O.OO78±O.OOO8M0yr"^Mpc"^) and 
twice that from radio observations. The SFR from local su- 
pernova surveys (0.0081 ^ q 0021 Afeyr'^Mpc"^) is consistent 
with the estimate from Ha. For the SFR from UV emissivity, 
our result means (Eb-v) < 0.48 with the standard extinc- 
tion law. The increase of SK statistics by a factor of 5 should 
positively detect the supernova relic neutrino events. 
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Figure 4. Neutrino spectrum expected at the detector after neu- 
trino oscillation. Uncertainties in the spectral shape are repre- 
sented by shading. 



The constraint derived here is about 5 times higher than the 
SFR obtained from Ha and far-infrared, and twice higher 
than the SFR from radio observations. 

If we would take the opposite case, i. e. the maximum 
allowed flux (90% confidence for the 1MB events) and nic — 
8Mq, the limit becomes stronger by a factor of 4; it nearly 
coincides with the SFR from Ha, and already overshoots the 
SFR from radio. 

We emphasize that large uncertainties in the model su- 
pernova neutrino fiux calculations are significantly reduced 
by empirical constraints derived from SN1987A. For instance 
an increase of rj is compensated by an increase of efi'ective 
neutrino energy or else by an increase of neutrino luminosity, 
so that the neutrino flux in the energy range that concerns 
us changes little. Neutrino oscillation generically enhances 
the high energy tail of the neutrino spectrum. Under the em- 
pirical constraint, however, this is absorbed into the change 
of other parameters. As a result the prediction of supernova 
relic neutrinos at SK is modified little. 

The lower limit of neutrino reaction rates at SK 

Assuming \og^|> = —2.15, which is the lower value of 
SFR from Ha allowed within the error range (or SFR from 
UV with Eb-v > 0.19) and taking our minimum neutrino 
flux estimate and rric = WMq , we obtain 

RsK > 0.40 yr"\ (12) 

A similar limit is derived if we take the lowest value of the 
supernova rate: 

i?SK > 0.44 yr"\ (13) 

The SK should see the supernova relic neutrino events if 
they increase the sensitivity by a factor of 5. 
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